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Abstract 

We show there is a non-null probability to produce neutrinos of the "wrong" type in general 
decays, within the SM augmented by the known three massive neutrinos and nontrivial mixing. 
Such effect is equivalent to an intrinsic flavor violation (lepton flavor violation) at creation without 
, requiring neutrino propagation. The exact amount of flavor violation depends on the neutrino 

cni : 

flavor state to be detected. For realistic conditions, the violation is tiny but much larger than 



other indirect lepton flavor violation processes not involving neutrinos as final states. For neutrinos 
produced in a continuous spectrum, we show the effect is larger, in relative terms, for neutrinos 
produced in the low-energy portion of the spectrum. For muon decay, if the usual neutrino flavor 
Qh' state is assumed a flavor violation as large as 1% is possible in the channel /x^ — t- Ure'^i'e- We also 



discuss the relation between flavor violation and flavor indefiniteness for neutrino flavor states. 
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I. INTRODUCTION 



Family lepton numbers or lepton flavors {Lg, L^, Lr) are conserved in the SM to a great 
extent. The clearest sign of nonconservation appears in the neutrino sector where neutrinos 
where observed to undergo flavor conversion as they propagate long distances. We have 
inferred from such phenomenon that neutrinos are massive and mix themselves through 
large mixing angles. As such, lepton flavors are not exactly conserved quantum numbers 
within the SM augmented by massive neutrinos jl- J]. Nevertheless, the amount of indirect 
lepton flavor violation induced by neutrino masses, such as in — )■ ej, is extremely small 
and unobservable [Si!] in comparison to the direct effect of neutrino oscillations. 

Structurally, the existence of replicated families of quarks and leptons by itself is an un- 
explained feature of the SM. As far as electromagnetic and strong interactions are concerned 
each particle with the same electric charge behaves in the same way. Very different mass 
scales (and then the interaction with the Higgs), however, distinguish the three known fam- 
ilies. The presence of weak interactions introduces the quark mixing responsible for flavor 
violating processes. Quark flavor, however, are identified through the different quark masses, 
despite the impossibility to observe them freely. For massless neutrinos, the definition of 
neutrino flavor is intimately connected to the definition of lepton flavor, which can be defined 
exactly as a conserved quantum number, at least at the classical level. With the presence 
of tiny neutrino masses, lepton flavors can still be considered as approximately conserved 
quantities and neutrino flavors correspond to the superpositions of mass eigenstates (fields) 
that carry the lepton fiavor in the massless neutrino limit. In that respect, the definition of 
neutrino fiavor is unique within the SM in the sense that it is defined as a superposition of 
mass eigenstates. Such structure is responsible for neutrino fiavor oscillations. 

We intend to consider here a slightly different type of neutrino flavor violation, i.e., an 
intrinsic flavor violation that could be present without propagation. As neutrino flavors 
are deflned, in some approximate way, as certain superpositions of massive neutrino states, 
some level of flavor indeflniteness is expected at least at the order of (Am/E^) . If a 

strict deflnition is considered we can explicitly calculate the amount of flavor violation that 
arises. Such task was undertaken previously for neutrinos created in pion decay 0. In that 
work we have shown that the probability for the flavor violating channel vr — )■ /iz/^, /3 7^ /i, 
is given by 

V{n^^^u,)^lsin^2e(^y. (1) 



2 \2E^r 

The largest effect was found for /3 = r for which the probability was of the order of 10^^ 
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Although very small, such effect is much larger than the branching ratio of indirect flavor 
violating processes such as /i — )■ 67. 

In this work, we want to generalize such treatment to calculate the amount of flavor 
violation in general decays, in special, in muon decay. The generalization occurs in the 
sense of considering decays involving three or more decay products, since pion decay is a 
two-body decay that was kinematically easier to treat. This generalization is important in 
two respects. Firstly, pion and muon decays are the major sources of terrestrial accelerator 
neutrino experiments and atmospheric neutrinos. Second, it is important to check how 
the neutrino flavor violation effect changes with a continuous emission spectrum for the 
neutrino. For that reason, the emission spectrum for the production of the wrong neutrino 



flavor is also ca^ 
experimentally 



culated. It is known that the emission spectrum for z/g in muon decay agree 
8] with the SM prediction [ol and the effects of neutrino masses introduces 
negligible distortions. 

Ultimately, this work intends to clarify the validity of the usual definition of neutrino 
flavor and its dependence with properties such as localization aspects that is known to play 
a crucial role in neutrino oscillations Il0| • Recent discussions on the oscillation of Moss- 
bauer neutrinos [11] or entangled neutrinos [l2j shows localization is an important aspect 
and should be carefully analyzed. We are concerned, however, with a related but more 
fundamental issue, i.e., the definition of neutrino flavor and its universality. 

The outline of the article is as follows: in SecHIlwe derive the general formula for neutrino 
flavor violation in decays. In Sec. lIIIt we calculate the neutrino flavor violation in muon 
decay. In Sec. lIVI we define equal-energy neutrino flavor states and calculate the flavor 
violation probability for them. Discussions and conclusions are presented in Sec.|Vl The 
appendices contain auxiliary material that were chosen to be separated from the main text. 



II. NEUTRINO FLAVOR VIOLATION IN DECAY 

Within the SM, ordinary neutrinos can be produced from two elementary processes: (a) 
charged current processes (involving W^) and neutral current processes (involving Z^). For 
the latter case, one neutrino and one antineutrino (two neutrinos) are produced in accordance 



to lepton number conservation in the SM but still neutrino oscillations might be possible 13|. 
We will be interested here only in neutrinos produced through decays induced by charged 
currents. 

We will assume through this and the next section that neutrino and antineutrino flavor 
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states are approximately well described by the superpositions 



m 



^ U^,\u,{k)). (2) 

where |i/j(k)) and |z/j(k)) are well defined asymptotic states with definite masses [l^ . nor- 
malized as (k|k') = 6^(k — k'). The latter normalization will be used for all momentum 
defined states henceforth. The orthogonality among the mass eigenstates |^'j(k)) implies 

Mk)\up{k')) = 6o^p6%k-k'). (3) 

Furthermore, we are interested in assessing the probability of lepton fiavor violation for 
neutrino fiavor states ([2]) produced in decays. We can distinguish two types: 

(A) l^^p^ + X, (4) 

(B) I^l^ + u^ + X, (5) 

where X is one or more particles without net lepton number. There is lepton fiavor noncon- 
servation, in the sense of an approximate family lepton number, if /3 ^ a. Types A and B 
of fiavor nonconservation are classified according to the fiavor nonconservation between (A) 
an initial and a final particle or (B) between two final particles. Pion decay is an example 
of type B decay without the accompanying product X. The neutrino fiavor violation in 
that context was calculated previously [3] . Muon decay is an example where both types 
of violation may occur and it will be treated in Sec. lIIII Usually, however, only one type 
of violation will be possible. The CPT conjugate processes of (jl]) and (E]) can be equally 
considered. 

Because of the definitions in Eq. ([2]), the amplitude for these decay channels in Eqs. (jlj and 
(|5]), with fixed a, should be regarded as a coherent sum of the amplitudes of the 3 channels 
involving the neutrino mass eigenstates Ui, i.e., (A) — )• z/j + X or (B) /—>/+ + z/j + X, 
i = 1,2, 3. The weight of each channel is dictated by the SM weak interactions, considering 
non-trivial mixing and non-degenerate masses for the three families of neutrinos z/,. 

Once we have the time evolution of the decaying state \la{t)) or \I{t)), the fiavor violation 
probability can be calculated by 

Vi^^u,{t) = l[d^Pp]\{x,up\Lm' = j[d^v,]\Y,up,{x,v,\Lm\ (6) 

'P^l.u,{t) = J[d'PF]\{X,L,l.p\im' = y"[rf^PF]|5^f//3,(X,Lz/,|/(t))P (7) 
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where Jld^Pp] denotes the integration over all final momenta. The sum over final spin states 
is implicit. We can associate a wave packet ipj to the parent particle, encoding informations 
on momentum distribution and spatial localization, by setting as the initial states 

(8) 
(9) 

We can assume that for a time t ^ l/F, where F is the decay width of the parent particle, the 
parent particle decays, among other possible states, into the final states containing neutrinos 



\la) = J d^piJl{p)\la{p)) , 

I J) = / d'pMp)\i{p)) ■ 



[[d'PF]Xaj{{PF})e-'^^-'\X,uy,{PF}) , 
J2 j [d'PF] Xaj{{PF}) e-^^-^*|X, L, u,; {P^}) 



(10) 
(11) 



where Epj is total energy of the final state, Pp denotes the final total 3-momentum and 
{Pf} is the set of the final momenta. 

A calculation based on the Wigner-Weisskopf approximation 0, [ist yields 



U*-F ■ 



N, 



-1/2 Kjf^aj 



iA 



(12) 



where we used the shorthands AE^j = Ej — Epj, N^j = {2u:)^'^2Ei[2Ef\ and A = ^, with 
r being the total decay width and 7 = the Lorentz factor; \2Ep\ = Y\^{2TtY2Epi is 
the product of the final energies; 1 = 1^ for type A decay. The factor AE^j — iA in the 
denominator is responsible for ensuring approximate energy conservation since 

A 



A-s>0 



lA^-zAP 



n5{AE) . 



(13) 



The factor A4aj is the invariant amplitude for the momentum defined process, with the 
mixing factor U*j factored out 16|], while / is a form factor necessary to regularize the 
expressions for large momenta 7]. To compute the fiavor violation probabilities, however, 
the form factor is not necessary [7|] . Recall that the production probability density for each 
mass definite channel is given by IXajP- 

Taking Eq. f lT2|) into account, the creation probability for both type A and B decays are 



J d'p mv)?j [d'-PF^iPF - p) I ^P^' 



U^- F ■ 



(14) 
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provided that we use the appropriate quantity Fr.i. We see the exponential e"*'^''^ * is respon- 
sible for the neutrino oscillation phenomenon 0, jo]- Despite of that, we can see the sum of 
probabilities 



f^=e,fi,T 



is time independent and will be shown to be approximately equal to the total probability 
of neutrino production. Ultimately, however, we will be interested in the flavor violation at 
creation, i.e., at times t that satisfy l/F ^ t ^ Lose- We will assume henceforth that the 
time in question satisfies such regime and references to time will be suppressed. 

To check the correct normalization explicitly, we consider the probability in Eq. fll4p for 
the type B decay / — )■ la^pX and take the limits: (a) massless neutrinos, nij — > 0, (b) small 
width f lT3|) and (c) parent particle at rest and with sufficiently small momentum uncertainty, 
i.e., |'?/'/(p)p is only appreciable around p ~ 0, within a size |p| < dp, for which the rest of 
the integrand varies very slowly. We obtain 



vii ^rx) = j ^^np, - p.) \M.\' = • - ' , (16) 

where P/ = (M/, 0) and Eq. flTB]) coincides with the branching ratio for a = /3. The nota- 
tion 1/'°' indicates we are considering the respective neutrino, or combination of neutrinos, 
massless. Notice Eq. (1161) is Lorentz invariant and considering the parent particle at rest in 
restriction (c) above is not essential. Notice, the order of the limits is important and we are 
taking (a) before (b); if we take (b) before (a) we obtain the incoherent limit discussed in 
Sec.|Vl It is also important to emphasize that Eq. (ITBl) is flavor diagonal (~ Sap), thus con- 
firming that neutrino fiavor is a well defined concept for massless neutrinos in the SM 
The same analysis can be performed for type A decays. 

For comparison we can calculate, under the same conditions but finite neutrino mass, the 
probability for the mass definite channel 



We can define the total probability 

Vtot{i ^ LvaX) = ^v{i^ LvjX) . (18) 

3 

In special, when there is only one channel involving one charged lepton such as e'^u^iy^, 
the probability above is unity. The probability (fTSjl is the usual production probability when 
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neutrinos are not detected 18|]. When the emission spectrum is considered for the charged 



lepton, the distortion of the endpoint compared to the massless neutrino case gives us the 



effective absolute neutrino masses 



18 



19| . In general, however, we can approximate 



Vtot {I ^ LvaX) ^ V{I ^ Ly^^^X) . (19) 

The difference is negligible [~ (Am/i?^)^] and goes to zero for massless neutrinos. Then, the 
sum of probabilities in Eq. f llSp is equal to Eq. f lTSjl within the small width approximation. 
If we divide Eq. (1151) by Eq. f lTSj) we obtain 

E = 1 ' (20) 

l3=e,ii,T 

where P;,^^ = Vi^up/{Vi^,yJtot- 

To proceed further in the analysis of Eq. (1141) . we can consider some approximations. Due 
to Eq. (1131) . approximate energy conservation holds and the amount of violation is of the 
order of F. Then, around the energy conserving values, due to simple kinematics, neutrino 
masses are negligible in the terms Aiaj and Naj which allows us to approximate them to 

Ma^{M^,)l%, N^^{N^,U^^o, (21) 

where EC denotes that energy conservation is strictly assumed |20i] . In other words, neutrinos 
are always produced ultra-relativistic and neutrino mass is not the leading term within any 
multiplicative factor. Approximations in Eq. (12 ip allows us to write Eq. (I14p as 

p..^/^i.,(p)r/gf(2.f.'(p.-p)i>..iiE^r. (-) 



Moreover, within the approximations of Eq. (I2ip . it can be shown (appendix [X]) that the 
integral in p can be decoupled and we can rewrite Eq. ([22]) as 

where all the quantities are calculated in the rest frame of the parent particle /. 

We can also make the flavor violating contributions explicit by rewriting the term inside 
the square modulus in Eq. ( IT^ as 

3 3 

UajU;^F^j = 6^^F^i + J2 UajU;^AF^j , (24) 
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where AFq,j = Faj — Fai- Thus the square modulus becomes 

3 3 3 2 

j=l j=2 j=2 

We recognize that only the last term of Eq. (125|) is flavor non-diagonal. If approximations 
are considered, we can write AFaj as 

AF ~ (^^-^^ 7 - 2 or 3 (26) 

where AE^,^ = E,y^ - . 

Specializing to a ^ P, under the approximation of one dominant contribution (two fam- 
ilies), the initial creation probability yields 

Vi^u, = \UajU;^\' j \Mp)\'j [d'PF] |AF„,f , J = 2 or 3. (27) 



III. NEUTRINO FLAVOR VIOLATION IN MUON DECAY 

We will calculate here the intrinsic neutrino flavor violation probability for muon decay 
/i"^ —7- z/ae"*"!/^. Both types A and B of violation will be possible, i.e., (A) a ^ ^ and (B) 
13 ^ e. The calculations developed here can be generally adapted to consider the flavor 
violation probability for any three body decay with neutrinos emerging as final states. 

In terms of mass eigenstates, the decay /i"^ — )■ Ua^^^-'is should be viewed as a co- 
herent superposition of the 6 channels /i^ — > ViC^Vj, i,j = 1,2,3. The calculation 
for each mass eigenstate channel with fixed momenta follows the momentum convention 
fi+{p) Ui{k)e+{q)uj{k'). 

Following the description of SeclITl the initial muon state 

|/i) = j rf3p^(p)|;,(p)), (28) 

after a time t ^ l/F, decays into the state 

J2 I rf'qrf'krf^k' x^.,ei(q, k, k') e-*^--*|e(q)z?i(k)z/,(k')) , (29) 

where Epij = Se(q) + ^p,(k) + E^^(k'). Considering Eq. (JED, X/.i,ej(q, k, k') = 
U^iU*jFfj^i^ej4'(PF) and the initial creation probability for yU"*" — > UaC^up is then 

r^9^>^, = j d'p mp)\'J d'qdVk'5'{PF - p) I E u,iU:,upju:^F,,^ejf , (so) 
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where Pir = q + k + k' and 



TJ TT* p 
'-^ fii'-^ ej-'^ 



-1/2 



AEij - iA 



= GU^,U:^M^j, 



(31) 
(32) 



where G = 2^20 f, = E^{p) - Ee(q) - E^^k) - E^^{k') and Mij is given in Eq. f lD3|l . 
The probabihty amphtude squared for unpolarized muon is 



l{p-k'){q-k) 



(33) 



spins 



By using the approximation Eq. (1501) becomes 

r[d^PB 



(34) 



where the muon is at rest. 

After following the calculations described in appendix [Bl we can obtain the flavor violation 
probabilities for type A and B. 

The type A neutrino flavor violation probability, within two families approximation (j27|l . 
is given by 



* |2 



dXjySp {xu)pii{x„) , 2 = 2 or 3, 



(35) 



where x^, = Ep/Wp, Wp ~ M^/2 is the maximum antineutrino energy (0 < X;^ < 1) and Sp^ 
is the emission probability (energy spectrum) for i/^ in the ordinary channel fi~^ — )■ v^e^Ve'-, 
see Eq. (IDlip . The flavor violating spectrum is then modified by the mixing factor in front 
of the integral in Eq. ( 15^ and 

{^Epf 



(AEpJ2 + 4A2 



which depends on Oi = nii/Wp, aij = |(aj + aj), 



(36) 
(37) 



1 1 
mf — m- 



2TW,, 



bij{x) 



x^ + af \/ x'^ + a'j — X' 



(38) 
(39) 



The monotonically increasing function bij{x) has range [0, (Aa)^/4], where Aa 
thus afj + bij < \{af + 4m^/M^ which is negligible. Such function is similar to the 

function / found in Refs. 2l|, |22| in another context. 
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The distribution function iSp^(x,y), at tree level and null electron mass, can be found in 
appendix It is a monotonically increasing function going from to a maximum. On 
the other hand, pij in Eq. fl36|) increases for small x. Within these approximations, the 
expression f l35|) can be integrated and we obtain, at leading order f lB3|) . 



P(/i+^z^„e+z.f ) = ^\U^iU*J^ul , 2 = 2 or 3. (40) 

We have typically 

V{^i+^yee+v'i'') = 8|f/^2f/e2rt^2i ~ O-Sw^^s ~ 5 X lO"*^ , (41) 

V{ji+ ^ P.e+i/f ') = 8|f/^3f/;3l'w|i ~ 2 culs ~ 10-2 , (42) 



where we have used |f/^2f^e2p ~ I^A»2t^r2p ~ 0.1, lU^sU*^]'^ ~ 0.24 [23| and we have assumed 

I^M3f^e*3lV|f^M2f^e*2l' < (A^?2) '/(A"^^)' ~ lO'^ in Eq. 1^ ■ If l^/^at/^gp ~ 10-3, we have 
to consider, instead of Eq. ( |4T]) . 

P(/i+^z^ee+i/i°') = SlUf^sU^sl^ojl^ ~ 0.8 x 10"^ ~ 5 x 10"^ . (43) 

We made use of the experimental values 

^2 f^^y^5 8^10-^ u;l,= (4^Y ^5.7xl0-\ (44) 

^2 v2iy^r^/ ' 23 v2iy^r^/ ' ^ ^ 



with lAm^gl 7.6 x lO^^eV^ and lAm^gl ^ 2.4 x 10~3gY2 



231 ]. Additionally, if we had used 



Eq. flB3p for UJ23, we would have obtained a value 18% larger in Eq. fH2|) . 
Analogously, the type B neutrino flavor violation probability is given by 



Jo 



(45) 



where Xi, = E^/W^, ~ M^/2, and the remaining functions are the same as in Eq. f l35l) : 
the energy distribution Su^{xl) is different and can be found in Eq. flD12l) . After integration 
we obtain, at leading order, 

Vif,+ ^i?ll\+up) = 12\U^,U:^\'u],, J = 2 or 3. (46) 

The typical values for /3 = /i and (3 = t are 

P(/i+^i^<"'e+i/^) = 12|f/^2f/:2l'w2i ~ 1-2 Wi'2 ~ 7 X 10-*^ , (47) 
Viii+^u'i^^e+ur) = l2\Ur2Ul2?ujl-^ ~ 1.2^^^ ~ 7 X 10"^ (48) 
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where we used |f/r2t^e2p ~ 0-1) addition to the values aheady used after Eg. (H^ . If 
|f/^3f/*3p ~ 10"^, we have to consider, instead of Eq. ( Hgj) . 

P(/i+^z^;°'e+i^^) = 12\Ur3U*^\'^ul^ ~ 1.2 X lO^^w^g ~ 7 X 10~^ . (49) 

In addition, if we had used Eq. (IB3p for CO23, we would have obtained a value 25% larger in 
Eq. (US]). 

The difference in the factor 3/2 between type A violation probability in Eq. (140 p and 
type B violation probability in Eq. fj46|) can be understood by analyzing the overlap between 
the function pij in Eq. fl36|) and the emission spectra of the would-be z/^ (type A), Su^{xu), 
and the would-be (type B), 5y^(x^). These distributions can be found in appendix [Dl 
Considering that Pij{x) is larger for small x, the overlap with S^^{x'j^) is larger than with 
Sp^X'^v) because the former is centered around a smaller value. 

It is also important to emphasize that from Eqs. (1350 and (H5|) we can extract the flavor 
violation emission spectrum for neutrinos in type A violation, /i+ — )■ Va^^^e-, a 7^ /i, or 
type B violation, /i"*" — )■ z/^e"^z/^, /3 7^ e. The emission spectrum of the wrong-type neutrinos 
relative to the spectrum of the correct-type neutrinos is given by the function p{x) in Eq. (1371) . 
Such function increases for smaller neutrino energies, which implies that the relative effect 
of neutrino flavor violation is larger for neutrinos produced in the low energy part of the 
spectrum. 

IV. EQUAL-ENERGY NEUTRINO FLAVOR STATES 

By assuming the usual neutrino flavor states ([2]), we have found in Sec. lIIII a flavor 
violation as large as 1% in the channel p^ — )■ VT-e^Vf. relative to the dominant channel 
/i+ — v^e^Vf.. It is then mandatory to check if a different definition of the neutrino flavor 
states can minimize the inherent flavor indefiniteness that arises. 

Qualitatively, we know from the relevant quantity uij (138|) that flavor violation was 
present because the size of the neutrino energy differences AEy ~ /S.m^/2Ei, is compa- 
rable to the intrinsic energy uncertainty of the creation process given by 5E ~ F. Thus, 
it is important to remark that although the calculations took into account both the energy 
uncertainty quantified by F and the momentum uncertainty o"p encoded in the parent par- 
ticle momentum wave function ^/'(p), the latter did not play any relevant role in the flavor 
violation probability (it could indeed taken to be zero). The reason lies on the fact that by 
calculating the probability to detect the neutrino flavor states ([2]), we are implicitly summing 
over neutrino mass eigenstates with the same momentum. 
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We can define, instead, the equal-energy neutrino flavor states 



|z/,:E,k) = f/,i|i?,(k,k)), k, = ^E^ - m\ . (50) 

The superpositions involve neutrino states with the same energy since E^-{kik) = E, i = 
1,2,3. It is also necessary to adopt the convention E > = max(mj). Notice, the states 
fl50|) do not obey the orthogonality condition ([3j). 

We can recalculate the flavor violation probability of Sec. fillip for muon decay by adopt- 
ing Eq. flHJ]) . For the type A violation yU+ — > UaGi'e^, by projecting the state fl2^ into 
(z/„:E,k|(z/(°'(k')|(e(q)|, Eq. is modifled to 

P(/i+^P,e+i.f ) = f d\{d^k)Ed'k' I J2 U^^rK^F^M^. + k.k + k')|' , (51) 

where / ((i^k)^; denotes dEE^ J dVLk- (Other choices are possible but they are not rele- 
vant as long as approximation fl2T]) is valid.) We have used 

(z^,:E,k|z7,(k')) = f/;,53(k'-k,k), k, = sjE^-mf. (52) 

We then use approximation (1^ to rewrite 

P(M+^^„e+.f ) = Jd'pd\{d'k),d'k' S\p - PfJ^I E ^l^r ' (53) 

where P^i = q + kik + k', Ak^ = k^ - ki, {d^k)i = ^^(kik) and 

ijji = i/jip + Akik), (54) 
AE^i = E^ip + Ak,k) - Ee(q) - E - E,ik') . (55) 

In the two families approximation, Eqs. (15^ and ( 155|) contribute as 



AE^i-zA A£;^i-iA A£;^i-iA VAE^^-iA 

A^i =^Pi-^|Jl ^ Aki k-^^(p) , (57) 

1 \ ^ ^^(P)-^^(P + Ak,k) 
UE^,-zAy (A^^i-2A)(AE^,-zA)- ^ ^ 

Let us now calculate the flavor violation probability (!53|) in two regimes: (i) ^ 
1/^r^M^ [lAV^il dominates over |E^(p) - E^(p + Ak,k)|] and (ii) < 
\\E^{p) — -E'^(p + Akjk)| dominates over \A^i\]. The details are shown in appendix [Bl 
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fi) a, > 



P(/i+^z^„e+z.^°') = l\U^iU*J\a,W,fj^dxS^^{x)yx^ + \Aa^\-xy (59) 



= W^^^*c.^\\-^^) , ^ = 2,3; (60) 

(ii) a, < i/v/ryw;: 

P(^+^P„e+.f)) = |f/,.f/:,p(^)'^rfx5p^(x)(/;;^7^-x)' (61) 

2 

= 4\U,,U:ful^ , 2 = 2,3. (62) 

We have used crl = J d^p\dijj{p)/dp\'^ and dp = / (i^p|pp|V'(p)P for the position and 
momentum uncertainties of the parent particle, respectively. (Assuming the mean position 
and momentum are null.) The shorthand Aaf = of — a\ was also used. 

We can compare the contributions of Eqs. (plj) and ([S2D by using cXp ~ l/2(Ta; and inserting 
typical values for a = r and i = 3: 

(i) P(^+^p,e+i^(«)) ~ 6x10-13 f^)', (63) 

Vlcm/ 

(ii) p(^+^p^e+4o)) ^ 6 X 10-29(^1^)'. (64) 



We conclude that both contributions are negligible unless is extremely large or small 
compared to 1cm. In fact, Eqs. fIBUl) and flB2]) have opposite behaviors as functions of a^', 
they are comparable when Ta,j. ~ (jp/M^ or 

~ ^y^^ ~ 1/^"^ ' (^^' ~ V^p)- (65) 

In that case both contributions of Eqs. 0631) and flMl) are of the order of 10"^^. 

It is important to remark that with the typical values of Eqs. (jHS]) and flM|) . approximation 
(|2T|) can not, in general, be applied to Eq. (15T]) to obtain Eq. ( 153|) . because we might be 
neglecting contributions of the order of 

Am| ^^_i7/10MeVA2 
2^2 



10-"^— J . (66) 



Therefore, to properly calculate the flavor violation probabilities for the states flSOj) . it is 
necessary to perform the full calculation of Eq. f l5ip . Nevertheless, the calculations of this 
section show that the flavor violation probability is negligible for the states fjSOjl . 
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V. DISCUSSIONS AND CONCLUSIONS 



We have calculated the flavor violating creation probability for neutrinos produced 
through decay, using the full QFT formalism within the Wigner-Weisskopf approximation, 
assuming that the usual neutrino flavor states ([2]) describe appropriately the neutrinos pro- 
duced in nature. The neutrino flavor violation occurs at creation and no propagation distance 
is necessary except for the lifetime of the parent particle. The calculations in this paper 
differ from previous QFT treatments (see Ref. 2J]) in the focus on the intrinsic flavor vio- 
lation, with the considerations of the decay width and the realistic interactions responsible 
for neutrino creation. We also managed to calculate the flavor violating neutrino emission 
spectrum. 

For the particular case of muon decay, the amount of flavor violation might be surprisingly 
large as 1% in the channel — t- UrC'^Ue compared to the ordinary channel fi^ — )■ v^e^Ve- 
We have found that the amount of flavor violation in such channel is much larger than other 
flavor violating channels because the v^- mixing is large and the relevant mass difference 
Am^3 ~ Am23 is relatively larger. The relevant quantity is given by lOij ~ I\m\^j2EyY in 
Eq. Compared to neutrinos produced in pion decay 7r+ — )• fi'^u,,,, where the amount 

of flavor violation of the same type was at most of the order of 10~^ [7], flavor violation in 
muon decay is much larger because the decay width is much smaller, while the remaining 
variables such as neutrino energy are similar. 

Considering the amount of neutrino flavor violation calculated in muon decay, we should 
consider the limits imposed by the CHORUS 2^ and NOMAD [26| experiments that was 
searching for the conversion z/^ — )■ 1/^ but excluded such channel at the level of 10^^. However, 
most of the z/^ beam comes from pion and kaon decay; the flavor violation probability for 
neutrino states ([2]) originating from pions is of the order of 10~^ 0] which is not constrained 
by such experiments and, for neutrinos originating from kaons, the flavor violation should be 
slightly smaller due to larger decay width. On the other hand, there is a 5.6% of contribution 
of z/^ coming from muon decay in CHORUS 271]. If the sensitivity to detect Ur is the same as 
for Ur, the probability estimated in Eq. fj42l) would be excluded. Moreover, if we extrapolate 
our results in Eqs. fj77|l and fHHl) to the production of z/g through the usual beta decay or the 
decay of long-lived heavier nucleus, we obtain u^^ ^ 1 which corresponds to the incoherent 
limit fl^7|) . Since no flavor violation is observed in neutrinos produced in nuclear reactors, 
we have to conclude that the neutrino flavor states ([2]) do not describe appropriately those 
neutrinos. Therefore, in contexts as common as muon decay or beta decay, the equal- 
momentum flavor states ([2]) are not appropriate. 
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It is also important to discuss an important limit: the incoherent limit when F — t- (or 
Uij — i- oo). Taking such limit in Eq. ( E^i or Eq. fHH]) we obtain 

V,^,^ = 2\U,,U*J' or V^eu, = 2\U0,U:^\' . (67) 

We easily recognize the expressions above as the incoherent limit if we resort to the two 
families approximation where we would have 2|?7^j?7*jp = | sin^ 26, where 9 is the associated 
two-family mixing angle. The results above (1671) are expected because in the incoherent limit 
the vanishing energy uncertainty would destroy the quantum coherence necessary to create 
the flavor states ([2]) and each mass eigenstate neutrino would be produced incoherently 10]. 
The same conclusion can be reached if we analyze the incoherent limit in fl22p : the square 
modulus of the sum in the last factor would be equivalent to the sum of the square moduli 
because the mixed terms would vanish due to the lack of overlap; hence the expression would 
be equal to Eq. f|T5|) . Moreover, the flavor violation calculated here is not negligible exactly 
because the energy difference among the different mass eigenstates that compose the muon 
neutrino, ^.E^,- ~ Am1-/2Ey, is comparable to the energy uncertainty imposed by the decay 
width 5E ~ r. 

For that reason, we define equal-energy neutrino flavor states (l50i) in Sec. (1IV|) and calcu- 
late, for type A muon decay, the probability do detect such states summed over all energies. 
In this case, the correct amount of flavor violation depends on the position (momentum) 
uncertainty of the parent particle (~ l/2(jp). As expected, if CTp ~ F the intrinsic flavor 
violation in Eq. (1601) is, except for numerical factors, identical to the flavor violation for 
equal-momentum states ([2]) calculated in Eq. fH2l) . However, fXp ~ F is equivalent to cr^ ~ r 
which is macroscopic for muons, i.e., ct = 659m. Position uncertainty should be smaller 
than the order of 1cm, or even much smaller (10~®cm) if the muon decays in a medium 28|. 
(For Mossbauer neutrinos, the uncertainty would be of the order of atomic size Even 
for ax ~ Im, we would obtain a flavor violation probability of the order of 10"*^ and thus no 
appreciable amount of flavor violation is expected. Although, it should be remarked that 
such flavor violation probabilities are much larger than the ones for indirect flavor violation 
processes such as — ?■ 67, with branching ratio ~ 10^^°. On the other hand, we can con- 
clude that the equal-energy states (ISUI) describe more accurately the neutrinos produced in 
muon decay than the equal- momentum states fl50l) . the description being more accurate for 
smaller decay widths or longer lifetimes of the parent particle, e.g., for the usual beta decay. 

Analogously, in the recent controversy concerning Mossbauer neutrinos, it was shown from 
a careful theoretical analysis 111] that neutrinos oscillate despite the tiny energy uncertainty. 
The reason is that the momentum uncertainty, which can not be as small as the energy 
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uncertainty, should be taken into account. In that case, equal-energy neutrino states also 
describe more accurately the neutrinos propagating from source to detector, enabling flavor 
oscillations. 

Two quantities control intrinsic neutrino flavor violation in neutrinos created through 
decay: the decay width ~ F (energy uncertainty) and the momentum uncertainty ~ ap of 
the parent particle. The former is intrinsic to the parent particle while the latter might differ 
depending on the process of creation, e.g., decay in vacuum as opposed to decay in a medium. 
From the calculations performed we can conclude that the decay width is relevant only for 
detecting equal-momentum states ([2]) (dp can be taken to be zero) while the momentum 
uncertainty is relevant when equal-energy states (l50l) are considered. (The decay width 
entered in Eq. (1^ but the suppression factor makes it usually negligible.) Reference |29 
derives exactly the quantity inside parenthesis in Eqs. (pi]) and (1^ to be smaller than unity 
for flavor oscillations to take place (called ACC and SFC conditions, Eqs. (4.2) and (4.7) 



in Ref. 



29 



respectively). We have shown here that such conditions are indeed necessary 
to avoid intrinsic flavor violation and ensure initial flavor definition for equal-energy states 
fISU]) . Interestingly, the condition ^ 1 fl55]l . obtained here for equal-momentum states 
([2]), is not obtained in Ref. 2^. 

Considering that any other choice for the neutrino fiavor states other than (j2]) or (l50|) still 
induces different contributions for the distinct mass eigenstates with differences comparable 
or larger than the contributions coming from F or cTp, we can extrapolate that there is a 
minimum amount of flavor violation of the order of 

> |f/«.t^|J'(^^) , cr, = max(F,ap), (68) 

for types A or B decays, as long as a* ^ Wi^ and the quantity inside parenthesis in Eq. fIBS]) 
is much smaller than unity. The exact amount of flavor violation would depend on the 
details of the neutrino state being detected. 

We should also discuss two aspects of the same phenomenon in neutrino creation: (a) 
flavor indeflniteness and (b) flavor violation. Both effects are related through the overall 
conservation of probability (120 p . i.e., a non-null probability to detect the "wrong" neutrino 
flavor implies that the probability to detect the "correct" flavor should be deflcient by the 
same amount when com par ed to the usual result. These effects are expected, at least, at the 
order of {Am/E^y |30| pointing toward the impossibility to deflne the neutrino flavor in 
an exact manner. Such result can be derived simply in flrst or second quantized formulations 
of flavor oscillations j6,|2l|. It is possible to avoid this intrinsic flavor violation [sol by deflning 
an inequivalent vacuum and different flavor states (3l|, but other problems appear (32|. In 
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this work and in Ref. j?!, however, we have shown through concrete calculations that intrinsic 
flavor violation probability can be much larger, of the order of (Am^/ (Q.E^a^^))'^ (1681) . 

The dependence of the neutrino flavor violation effect on the decay width and momentum 
uncertainty brings about another possible effect: (c) the source dependence of neutrino fla- 
vor {^,[331 . The difference between the probabilities calculated through the equal-momentum 
states ([2]) and equal-energy states fl5U]) also indicates that neutrino flavor could depend on the 



detection process [29|, l33|. We did not pursue such effect here, choosing to focus on idealized 



measurements (summation over momenta or energies) which insured automatic probability 
normalization. In any case, the calculations presented here suggest that neutrino flavor is 
not universally defined and the effects might not be negligible as usually assumed [5i] . For 
example, for the same observable, i.e., the probability to detect the state |z/^) as defined 
in Eq. summed over momenta, yields the expected result at the level 10~^ in the de- 
cay 7r+ — 7- /i'^z/^ while the same deviates from the usual expectation by 1% in the decay 
fi^ — 7- v^e~i'e- Even with a different definition of the z/^ flavor state, large flavor violation 
would remain if we could hypothetically produce free muons (at rest) with large position 
uncertainties of the order of its lifetime r (or dp ~ F). 



One aspect that was not considered here was the role played by entanglement 121]. How- 
ever, even in that case, some level of flavor indefiniteness should be present since the mo- 
menta and energies of the neutrino eigenstates are determined to be distinct by the conser- 
vation of energy-momentum. Source dependence should also occur since a neutrino flavor 
state produced in one process might differ from another process by the values determined 
by the conservation of energy momentum. 

Incidentall y, s ource dependence as described here might account for the anomalies found 



in the LSND 



34| and MiniBoone 35| experiments. The flavor violation probability of less 



than 1% is the amount of violation necessary to explain the LSND anomaly. Although the 
required type of flavor violation in LSND (z/^ — )■ z/g) do not match the largest probability 
found in this paper {v^ — )■ z/^), some kind of source dependence could account for or, at 
least, ameliorate such anomalies without requiring any new physics beyond the known three 



neutrino families 
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, |37| . Moreover, the flavor violation effect found here agrees qualitatively 
with the anomaly found in MiniBoone (in neutrino mode) since, from Eqs. (!35|) . fH5|) . f l59|) 
and fl^T]) . the effect is (relatively) larger for neutrinos produced at the low energy portion of 
the muon decay spectrum. 

In addition, as a buy-product, if the amount of neutrino flavor violation at creation 
estimated in Eqs. fHl]) .f l42|) .f l47|) and fHHj) were detectable, we would gain an observable which 
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is extremely sensitive to ^13, i.e., V{fj,~^ — > Uee'^Ue). Depending on how large is the value of 
|?7e3| the contribution of Am"^^ might dominate over Amf^. Moreover, if |t/^3t/*3| ~ 10"^ 
the contributions from the two mass differences might be comparable and even CP violation 
might be observable. Of course, to probe such quantities it is necessary to detect z/g in the 
states fl50|) with precision better than 10^^ compared to the main channel. 

In conclusion, there could be neutrino flavor violation at creation for the muon neutrino 
produced through muon decay with detectable probability if equal-momentum flavor states 
(EJ are detectable in some way. In that case, the muon neutrino produced from pion decay 
might be slightly distinct of the muon neutrino produced in muon decay; the distinction 
being possibly observable. In general, however, equal-energy states ( 150|) describe more ap- 
propriately the neutrinos produced through muon decay and the decay of other long-lived 
nucleus, and no detectable neutrino flavor violation is expected, unless the parent particle 
has uncommonly large position uncertainties. 

Acknowledgment s 

This work is partially supported by the Brazilian agencies FAPESP and CNPq through 
grants 09/11309-7 and 309455/2009-0. The author thanks O. L. G. Peres for helpful discus- 
sions. 



Appendix A: Decoupling of Eq. ([22]) 

We will show here that Eq. f l22p within the approximation (I2T]) can be decoupled into a 
product of two factors, one of them being the decoupled integral J d^p\iljj{p) \ which is unity. 
In other words, we will show here that the factor 

does not depend on p. 

Firstly, the dependence on p of Eq. (lAip lies on ), |A^Q,p and A = r/27/, where 
7/ = Ej/Mj. Then, apply a change of variables {Pf} — ^ {Pp = ^J^Pf} such that the two 
sets of variables are related by a Lorentz boost A/ defined by {Ej, p) = A/(M/, 0). We can 
write 

^^P^ - 2E,(P) J 7/ ^ AE^, - ^T/2 ' ^^^^ 
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where the primed variables refer to the rest frame of /. We have also used the relations 
5^{P'p) = 77(p)(53(Pi7 - p) (notice p' = 0) and AE^ = -fiAE^j. Therefore, 

i.e., g{p) can be calculated assuming the parent particle / is at rest. We then obtain Eq. f l23l) . 



Appendix B: Some calculations 

We describe here some calculations necessary to get from Eq. fl34|) to Eqs. fl35l) and fj^Sjl . 
The details to obtain Eqs. (1591) and flM]l are also shown. 

For type A violation, Eq. (1351) . considering U/s massless and the two families approximation 
(127|) . we can rewrite Eq. (IMj) as 



r^,^.f^ = s.,\u,.u:f^^^^^ (Bi) 

where AE,^ = E.^^) - E,,{k), Ta"^ is defined in Eq. ([US]), Q = (M^ - ^p(k),-k), 
Ep(k) = ^{Ep-{k) + _E'p^(k)) and e = AEp^; Me is the electron mass. We can compare 
Eq. f lBip to Eq. (]D6|) . It is then possible to use the approximate expression (1C24P because 
the approximation is only inadequate for |k| ^ which corresponds to the endpoint 

of the neutrino spectrum. Near the endpoint, however, the emission probability is sup- 
pressed by the function pn. After an angular integration, a change of variables and use of 
Eqs. fIDllI) we obtain Eq. fl35l) . The function pij in Eq. fl36l) is obtained after the manipulation 
AEp, = (m^ - ml)/2E^. 

For type B violation, Eq. fH5|) . considering z/^ massless and the two families approximation 
fl27|) . we can rewrite Eq. as 

RC"^ f (AE 

where AE,^ = E,^{k') - E,,{k'), Q = {M^ - E,{k'),~k'), E,{k') = |(^.,.(k') + E,,{k')) 
and e = AEy.. Comparing Eq. (lB2p to Eq. (]D7p . after an angular integration, a change of 
variables and use of Eqs. (]D12p and (lC24p . we obtain Eq. (H^ . 
For Eqs. fj^Ol) and fj46l) it is necessary to use 

/ dx— = cj^fl — Va^ + g;^ arctan( , =)] 
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dx 



1 



2^2^-2 = X[l-(a' + a;')ln(l+ ^ 
oj^ + x^ + ^ + 



)]• 



(B3) 



To obtain Eq. flM?]) we firstly insert Eq. f l57p into Eq. f l53p and use the two families ap- 
proximation. Then, perform the integrals in k' and q by using T?^ in Eq. flC26p . We can 



assume any dependence of the integrand on p can be approximated to the central value 
p except for 



k- 



av'(p) 



dp 



1^2 

2^^ 



Equation is then identical to Eq. (ID6p with 
in the integrand, of the term 



A 



(B4) 

instead of T"'^ and the inclusion. 



(B5) 



The desired result is obtained after the approximations T^^ ~ T"^ and kj ^ ki followed by 
the angular integration in k and the integral in ki instead of E. 

Equation f l6T|) is obtained in a similar way. Firstly, insert Eq. fl58|l into Eq. fl53|l and use 
the two families approximation. Then, we perform the integrals in k' and q by using 7a"^ 
in Eq. f lCSp and take the central value p ~ except for 

kp 



-Hp) 



1 a: 



2 

p_ 



(B6) 



The additional term in the integrand is now 



1 (AA;0 



2 M2 2A2 



(TpW,\^ 1 



(mj 



m 



(B7) 



where the factor 2A^ comes from the definition of 7a "'^ when compared to T"^. After the 
approximation 7a ~ T"^ and the integration in ki we obtain the desired result. 



Appendix C: Two particle phase space 

The contribution of any two-particle phase space to muon decay can be quantified through 

T'^^{Q;m,,m2)^ J ^^^'^'(Q " " , (CI) 

where Ei = A/k^ + mf, = Ei, i = 1,2. The integral in Eq. (ICip can be carried out 
straightforwardly by following some steps: (1) rewrite = d'^ki5{kf — mf)9{k^), with 
free fcf, and perform the change of variables ki = A/2 + x, k2 = A/2 — x. (2) Exploit the 
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Lorentz covariance of the integral which ensures the property T"'^(A-^Q) = A"^,A'^.,T"''^'(Q) 



for a general Lorentz transformation A. (3) From the property in (2), write T°^((5) = 
+ /2(Q2)Q«g/3 and extract /i,/2 by calculating and T'^^QaQ^- The explicit 
form of T"'^ is found to be 



(C2) 



where M = mi + m2, Am = mi — m2 and 



(02 



(Am) 



2 1 



Af 2 (Am) 



Q2 



For massless particles, T"^ reduces to [38 1 



We can define a similar quantity for finite decay width, necessary in Sec.lD| 

2A3 fd%d% 6^{Q-ki-k2) k'^k^ 



TC 



(C3) 



(C4) 



(C5) 



2^x2^2 \Qo + ^e-Ei-E2-iA\'^\Qo-^e-Ei-E2-iA\^ ' 

where Ei = ^/k^ + mf, Q = {Qq, Q) is a four-vector, e is a number and the factors in front 
of Eq. (ICSp are so chosen because 

2 A3 



-> S{x) , as A 0. 



(C6) 



7r[x2 + A2]2 
In that way, 

7;"'^(g,0;mi,m2) ^T"^(Q;mi,m2), as A ^ 0. (C7) 

The calculation of the integral (ICSP is more involved than fICip because Lorentz covariance 
is lost but the step (1) described previously can be carried out and leads to 
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2A3 



4|Q| 



dxodXodip 



R 



IQo + l^-Ao-^AnQo-l^-Ao-^Ap 



(C8) 



where A" = (Aq, Q), x = a;||Q + x_l, is the angle in the x_|_ plane, 

x\\ = x-Q = |^(xoAo - /i^) , 
|x_l|^ = w^(Ao,xo, IQI) , 



u;^(Ao,xo,|Q|) = -—{xo-aY + ^yA^y), 



Q 



A2 



An. 



(C9) 
(CIO) 
(Cll) 

(C12) 
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The function A{X^) is the same as (lC3p . The region of integration R = Ri (1 R2 itl the 
(Ao, Xq) plane is constrained by 



i?i : co^> 0, R2 : 2|xo| < Aq . 



(C13) 



The constraint R2 comes from the ^-functions in step (1). We can also consider the equivalent 
constraints coming from the original variables 



i?3 : A^ = {ki + k2f > M^ Ao > , R^: {2xf = {ki - k2f < {Amf 



(CM) 



where M,Am are defined after Eq. flC2p . The constraints flC14p follows from the fact that 
the variables ki, k2 are forward time-like 4- vectors. 

Analyzing the constraints in Eqs. f lClSp and flC14p we conclude that the region of inte- 
gration involves 

Ao e [(Ac) 



where 



(Ao) 
(a^o) 



max 
min 



+ Q2 = Em(Q) 

a±i|Q|A(A2). 



(C16) 
(C17) 



It can be explicitly checked that the antisymmetric part is null 

' dxof l{x"X^ - X'^x^) = . 

After some lengthy calculations we obtain 

[{g'^'^X^ - A"A'^)A2(A2) + 3A"A^5(A2) 



(C18) 



rA"^(Q,£;mi,m2) = ^ / dXo A{X^^ 



A3 roo 
12 



(C19) 



'(Ao).„.„ |Qo + i£-Ao-^A^Qo-i£-Ao-^AP 

where A^ = Aq — and A,B are the same functions in Eq. (]C3[) . Since the integrand 
in Eq. (IC19P is only appreciable when Ao ~ Qo > (Ao)mm, provided that \e\ <^ A, we can 
approximate 

^Ao^go, if Qo >£;a/(Q), 

Xo = Em{Q), if Qo <^a/(Q), 

in all the terms except in the denominator and notice Em{Q) = (Ao)min- The remaining 
denominator can be integrated exactly and we obtain 

A2 



Ap(d): 



(C20) 



TA°''^{Q,e;mi,m2) 



A2 + l£2 



T^Hi^e, A, Qo-Em{Q)) T"^(Q; mi, ms) 



(C21) 
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where T"'^ is understood as the expression in Eq. (lC2p without the ^-function. The role of 
the ^-function is played by the function H which is given by 

1.1 , ^v-a\ . 1 



H{a, ^5 ?/) = 2 ~'~ 27r ^^^tan^^^j + arctan 



A 

Ana 



In 



(j/+a)2+A2 



(C22) 



Such function has the property that it is symmetric and localized around in a-direction 
and it behaves as a smooth ^-function (step function) in the y-direction. More specifically, 
H{a, A, y) is negligible if ?/ ± a < -A, H{a, A, ?/) ~ 1 - if y ± a > A, and H{a, A,y) ^1 /2 
if y"^ - <^ -A. 

Due to the properties of H, we can use the approximation 



Ap(e): ^' ^0 - ^^^(Q)) ^ A^i^^^^' ~ ' (^22) 



and obtain 



rA''^iQ,e;mum2) ^ ..^ , ^ T^^Q-mi^m^), (C24) 

A^ + 



where now the ^-function is included in T"'^ (IC2|) . The approximation Ap(e) in Eq. f lC23|) is 
inadequate only around e ~ and Qq ~ -E'm(Q) in a region of size ~ 2A. Notice Eq. (IC24p 
satisfies 7^°''(Q, 0; mi, ma) = T°^(Q; mi, ma). 

We can check the approximate expression in Eq. flC24l) is valid for Qq — Em{Q) ^ A by 
noting that we can extend the lower integration limit of Eq. flC19p to — cxd without changing 
the integral appreciably. 

We made use of the integral 

00 ^ 

dx T, tt; —rrr, r- 77^ = — TTT-f^fo, A, w) . (C25) 

_y [(x-a)2 + A2][(x + a)2 + A2] 2Aa2 + A2 ^ ' '^^ ^ ^ 

The integral can be performed by splitting the integrand into two terms, each containing 
exclusively one of the factors of the denominator. The integral in the limit y oo can be 
calculated explicitly by residues. 

We see the tensorial form of T^'^ (]C2ip is the same as T°'^ (\C2h within the approximation 

Ap(d) dmoD. 

For completeness, we also calculate the tensor 

Ta (g,mi,m,) = -y^^^^-^-^--^, (C26) 

where Ei = -\/kfT^7T|, Q = {Qo, Q) is a four-vector, and the factors in front of Eq. flC26p are 
so chosen because of the limit ([T5|l . The detailed calculation is similar to the calculation of 
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Th''^{Q, €] nil, rn2) fICSp above. The final result within the approximation Ap(d) in Eq. f lC20p 



is 

rA°^(Q; mi, m^) ^ T^^{Q- m,,m,)Ho{ ^''^l'^^^ ) , (C27) 
where T°^(Q;mi,m2) should be understood without the function 6 which is replaced by 

Hq{x) ~ ^ + ^ arctan(a;) . (C28) 



Appendix D: Muon decay 

Let us consider the muon decay n^ip) — )■ v ^{k)e^ {q)i'eik' ) , where the 4-momenta of each 
particle is explicitly written. 

Muon decay is described by the four-point Fermi interaction 

^ = -2V2GF Yl {la{x)rLU^,u,{x)Yu,{x)U]p^^Llp{x)) , (Dl) 

where L = |(1 — 75) and {Uai\ denotes the PMNS matrix. The invariant amplitude at tree 
level is given by 

M,, = 2V2GfU^,U:^M,, , (D2) 
Mij = D^(p)7"ivt;,,(k)n,^(k')7a^^e(q), (D3) 

while the square modulus, averaged over initial spin states and summed over final spin states, 
is given by 

spins 

The decay rate at rest is given by 

8G^ fd^q d^k d^k' 



T{fi^i?ll\i,'^') = -rl^-^ f ^^^6\p - Pf) ip-k')iq-k) , (D5) 

-^PahT"^{p-k-Me,0), (D6) 



2M^(27r)5 

8G2 r d^k 



r d-^k' 

J ^ip-k')n{p-k'-M^,0), (D7) 



2M^(27r)5 J 2E, 
MfGl 

where p = (M^, 0) and Pf = (g)e + ik)p + (A;%^. The tensor T"^ is defined in Eq. fICT]) . 
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We can rewrite Eqs. (1D6P and (ID 70 in terms of the energy distributions of and i^e, 
respectively, as 

/ dxp Sp (xp ) = Tf, dx^^S^^x^J. (D9) 



Jo Jo 
The energy spectra for e"*", i^e, t'/i at tree level, for rrii, = Me = are 



5e(Xe) = 2X^(3 - 2Xe) , (DIO) 

5p^(a;pJ = 2x^^(3 -2xpJ, (Dll) 
5,,(x,J = 12x^^(1 -x,J, (D12) 

where Xe = EjWe, x,^ = E,JW,^, Xp^ = E^jWp^, We ^ W,^ ^ Wp^ ^ ) is the 

maximum energy. The distributions (lD10p -( lD12p can be obtained from Eqs. (lD5p - (lD7p . 

Appendix E: Flavor diagonal term 

We intend here to analyze the second term in the r.h.s of Eq. f l25p which is flavor diagonal. 
We have estimated in Ref.0 that such term would be negligible for pion decay. However, 
Eq. (12 op indicates that such term has to be of the order of the flavor violating terms but 
negative in sign. 

For concreteness, let us consider the muon decay /i+ — )■ z/ae"*"//^"', neglecting the mass of 
one of the neutrinos. Let us also disregard all the channels except the dominant one involving 
the positron and neutrinos such that V = V in Eq. fl^Up . We then have for Eq. f l20p . 

V^^,^+V^^,^+V^^,^ = l. (El) 

Therefore, 

P^^p^ = 1 - V^,.,^ - V^^,^ , (E2) 

and the flavor conserving probability deviates from unity (Ptot) by the flavor violating prob- 
abilities calculated in Sec. lIIII 

For comparison, let us rewrite the second term in the r.h.s of Eq. fl2^ as 

2(AEpJAE,i 



"1 - 1^^^^ _ ^^|2|AE,i - ZA|2 

-{AE,y + 2AE,XE,-EF) ^^^^ 



|AE^i-iA|2|AE„i-iA|2 ' 

where AEp. = Ep- — E^^ and Ep = Ee + E (o) + Ep. is the mean final energy considering 
Ep^ = \{Ep- + -EpJ. Wee see the first term in the numerator of Eq. flE3p gives rise to the 
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same contribution, with opposite sign and different mixing matrix contribution, as Eq. fl27p 
considering Eq. (126 p . On the other hand, the second term in the numerator of Eq. (1E3P was 
estimated in Ref.lll (appendix E) and it is of the order of ufiiV / E^Y , which is neghgible. 
The same conclusion can be reached by making the exact calculation considering Eq. (ICSp 
with an additional Aq in the numerator. 

For muon decay, within the two families approximation, we can calculate the contributions 
from Eq. fl2^ explicitly, 

V{^i+^v^e+ur) = 1 - 8|f/^,|2u;2^ + 8|t/^,| Vi (E4) 
= l-8|f/^,[//3i|Vi, /3 = eorr. (E5) 

The three terms of Eq. (lE4p correspond to the three terms of Eq. (I25p . in the same order, 
and we see Eq. (lE2p is satisfied. 
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